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: Abstract 



This paper presents a complete set of quasilocal densities which describe the 



> 

\^ ' stress-energy-momentum content of the gravitational field and which are built 

■ with Ashtekar variables. The densities are defined on a two-surface B which 
O ■ 

Tjj- . bounds a generic spacelike hypersurface S of spacetime. The method used to 

! derive the set of quasilocal densities is a Hamilton- Jacobi analysis of a suitable 

I covariant action principle for the Ashtekar variables. As such, the theory 

' presented here is an Ashtekar- variable reformulation of the metric theory of 

O^' quasilocal stress-energy-momentum originally due to Brown and York. This 

^ . work also investigates how the quasilocal densities behave under generalized 

boosts, i. e. switches of the S slice spanning B. It is shown that under such 
boosts the densities behave in a manner which is similar to the simple boost 
^ ' law for energy-momentum four-vectors in special relativity. The developed 

■ formalism is used to obtain a collection of two-surface or boost invariants. 
With these invariants, one may "build" several different mass definitions in 
general relativity, such as the Hawking expression. Also discussed in detail in 
this paper is the canonical action principle as applied to bounded spacetime 
regions with "sharp corners." 
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I. INTRODUCTION 



The geometric expression for the energy of a nonrelativistic system (the functional form 
of the Hamiltonian in terms of the coordinates and momenta) can be discerned from the 
system's action functional This follows from a basic tenet of Hamilton- Jacobi theory: the 
classical energy of the system is minus the rate of change of the classical action (the Hamilton- 
Jacobi principal function) with respect to a unit stretch in the absolute Newtonian time. The 
ability to define the classical energy in this way rests on the fact that in the conventional 
variational principle for the system the lapse of absolute time is fixed as boundary data. 
From a practical standpoint, this means that in order to find the geometric expression for 
the system's Hamiltonian one need only consider the general variation of the action in which 
the endpoints of trajectories in the variational set are not held fixed (known as the Weiss 
action principle U|). Upon inspection of the boundary-term contributions to the variation, 
one can determine the canonical momenta as the factors which multiply the variations in 
the endpoint values of the coordinates. Furthermore, after the momenta are determined, 
careful inspection of the boundary-term factor with multiplies the variation in the absolute 
time then reveals the functional form of the Hamiltonian. 

Brown and York have proposed a generalization of the Hamilton- Jacobi method, which 
is applicable to a wide class of generally covariant field theories of a spacetime metric (in 
any dimension); and they have used this generalized method to discern what geometric 
expressions play the role of quasilocal stress, energy, and momentum in general relativity. 
Field theories of a spacetime metric enjoy a crucial feature in common with simple 
nonrelativistic systems: in the action principle it is possible to fix the time as boundary data. 
To see that this is indeed the case, consider a spacetime region which is topologically 
the Cartesian product of Riemannian three-manifold S and a closed connected segment 
of the real line /. The three-manifold S has a boundary dH = B (which need not be 
connected). Therefore, one element of the boundary DM. of is a three-dimensional 
timelike hypersurface T ("unbarred" T is reserved for a more special meaning) which has 
the topology oi I x B and is a (2 -|- l)-dimensional spacetime in its own right. The other 
boundary elements are t', the three-manifold corresponding to the initial point of /, and 
t", the three- manifold corresponding to the final point of /.[| Such a spacetime region 
is depicted in Figure (|I|). Now suppose that we are given a "suitable" action functional 
for the metric (and possibly matter) fields on the spacetime region Ai. By "suitable" 
we mean that the variational principle associated with the action features fixation of the 
induced metric on each of the boundary elements t', t", and T. In particular, the lapse 
of proper time between the initial and final hypersurfaces is fixed as boundary data since 
this information is encoded in the fixed T three-metric. The quasilocal energy is then 
identified as minus the rate of change of the classical action with respect to a unit stretch 
in the proper time separation between t' and t" . (Therefore, inspection of the boundary- 



^One may imagine that A4 C U, where U is some ambient spacetime known as the universe or 
sometimes the heat bath. The boundary B and its history T are simply collections of points in U 
and need not be physical barriers. 
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term contributions to the variation of the action can reveal the geometric expression for 
the quasilocal energy. This geometric expression is obtained by isolating the factor which 
multiplies the variation in the lapse function which controls the proper time separation 
between B slices of T.) However, notice that the T three-metric provides more than just 
the lapse of proper time between the initial and final slices, since it contains information 
about all possible spacetime intervals on T. One is free to consider the rate of change in 
the classical action which corresponds to arbitrary variations in the T boundary data. A 
quasilocal surface stress-energy-momentum tensor corresponds to this freedom. For the most 
relevant case of general relativity, the analysis of Ref. ||^ has demonstrated how this tensor 
leads to quasilocal surface densities for energy, tangential momentum, and spatial stress (all 
are pointwise tensors defined on B) which describe the stress-energy-momentum content of 
the S matter and gravitational fields contained within B. The theory of quasilocal stress- 
energy-momentum originally proposed in Ref. |@] is currently being extended considerably. 
One extension has been the introduction of quasilocal surface densities for normal momentum 
and temporal stress. The new developments associated with this extended theory will appear 
in an upcoming paper 0, and the results of the present paper are based heavily on these 
new developments (though the analysis here is reasonably well-contained). For a description 
of the new developments to be found in Ref. and how they relate to the present paper, 
see the discussion section at the end of this work. 

This paper uses a Hamilton- Jacobi-type method to derive quasilocal stress-energy- 
momentum surface densities which are built with the Ashtekar gravitational variables [Q, 
thereby fully extending some preliminary results concerning Ashtekar variables and quasilo- 
cal energy- momentum given in 0. Since the the Ashtekar version of general relativity is 
inherently a non-metric formalism, the Hamilton- Jacobi analysis given by Brown and York 
has to be modified. However, though the choice of Ashtekar variables over Arnowitt-Deser- 
Misner (ADM) variables @] leads to non-trivial technical complications and some conceptual 
ones, our overall goal is the same as in the metric theory; and the cornerstone of the method 
used here remains a "suitable" action principle, i. e. information about the lapse of proper 
time must be fixed as boundary data. Now, the usual covariant formulation of the Ashtekar 
variables is based on the well-known chiral action independently given by Samuel [0] and 
Jacobson and Smolin |^. This is a Palatini action which features the independent variation 
of the spacetime self-dual spin connection and the SL{2, C) soldering form. Applied to our 
spacetime region M., this action principle does not feature fixation of metric data on T, and 
hence it is not well-suited for our purposes. Perhaps, one could consider adding the necessary 
boundary terms to the chiral action in order to obtain a suitable variational principle. How- 
ever, here we follow another route which is based on a lesser-known covariant formulation of 
the Ashtekar theory which has been given by Goldberg. Goldberg's action functional is 
first-order, but in the variational principle the connection is not varied independently from 
the tetrad. We find that, subject to certain gauge fixation of the tetrad, Goldberg's action is 
a tetrad version of the action functional used to derive quasilocal stress-energy-momentum 
in the metric scenario. It should be mentioned now that partial gauge fixation of the tetrad 
and triad plays a crucial role in what follows. At first sight this may seem objectionable. 
But one should recall that such gauge fixation is also unavoidable in the triad formulation of 
Hamiltonian gravity when one discusses the notions of total energy and momentum in the 
asymptotically-flat scenario. In that case one must deal with a "fiducial triad at infinity." 
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0^ The gauge fixation of tlie triad in tlie quasilocal context is of the same nature. 

There is a subtle interpret at ional issue concerning the analysis to follow which deserves 
so comment at the outset. The Brown- York quasilocal densities are not unique, since one has 
the freedom to add a subtraction term —S^ (a functional of the fixed boundary data) to the 
gravitational action 5*^ which is used to derive the densities. Brown and York have offered 
the interpretation that such freedom allows one to set the reference points for the quasilocal 
densities. Now, the results of gravitational thermodynamics are, in fact, independent of the 
choice of subtraction term, and, therefore, such freedom seems to be an unnecessary one when 
examining the statistical mechanics of the strong gravitational field. [|r3,|T5] However, the 
subtraction term plays an important role in several other theoretical contexts. For instance, 
it must be incorporated into the definition of the quasilocal energy, if in the suitable limit the 
definition is to agree with the ADM notion of energy at spacelike infinity. [^,0 Furthermore, 
recent research has suggested that there is an implicit reference point set in some spinor 
constructions of gravitational energy based on the Witten-Nester integral. In this paper 
the passage from the triad ADM variables to the Ashtekar variables is effected by the addition 
of a purely imaginary boundary term to the action. We formally treat this boundary term 
as a subtraction term —S^ a la Brown and York. This allows us to construct the theory in a 
parallel fashion with the presentations given in Refs. However, though technically this 

viewpoint is completely satisfactory, it should be realized that it is less satisfactory from an 
interpretational standpoint. Indeed, if we wish to adopt the Brown- York interpretation for 
the imaginary subtraction term, then we are confronted with the issue of imaginary reference 
points for the quasilocal densities. Furthermore, even with the imaginary subtraction term, 
in the suitable limit the Ashtekar-variable expression for the quasilocal energy as given here 
does not agree with the ADM notion of energy at spacelike infinity. This seems alarming, 
but in fact is not a real problem. It merely signifies that it is perhaps better to view the 
imaginary subtraction term not as a true subtraction term, but rather as part of a different 
base action S'^ = — suitable for the Ashtekar variables. To derive an expression 
for the quasilocal energy in terms of the Ashtekar variables which is in agreement with the 
ADM expression, we would need to begin with an action which differs from this different 
choice for the base action by yet another subtraction term —5'"' (so the full action would be 
5" = 5"^ — 5""). In other words, our analysis is actually performed only on a base action 
S'^ (even though we split this base action into two pieces and treat one piece formally as a 
Brown- York subtraction term), and it should be understood that in some contexts it may 
be necessary to consider the addition of appropriate subtraction terms to this base action. 
We discuss these issues in more detail in the concluding section. 

The organization of this paper is as follows. In § II, the preliminary section, we discuss in 
detail the geometry of A4 in terms of several classes of spacetime foliations. This discussion 
is the groundwork for the analysis in the main sections. We also collect some notations and 
conventions in this section. In § III we derive a full set of quasilocal densities which are 



^This is quite analogous to the situation in nonrelativistic mechanics, where one can affect the 
definition of a system's energy and canonical momenta by adding boundary terms to a system's 
action. 
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expressed in terms of the Sen connection and triad on S, and thus may be easily rewritten 
later on in terms of the canonical Ashtekar variables. The geometric forms of these densities 
are discerned from a careful analysis of the boundary terms which appear in the Goldberg 
action principle. This analysis is quite analogous to the method described for nonrelativistic 
systems in the introductory paragraph. In § IV we turn to the issue of how the collection 
of quasilocal densities behave under generalized boosts. This behavior is similar to the 
simple boost law for energy-momentum four-vectors in special relativity. With the derived 
boost relations we then show how to obtain a number of two-surface or boost invariants, 
one of which is the Hawking mass.Q [131 Also in § IV, we consider the canonical form of the 
action principle for spacetime regions with "sharp corners." This analysis supplements recent 
results from standard metric gravity for such spacetimes. The appendices provide some 

kinematical results necessary for the central discussions. The first three appendices develop 
the results necessary to write down the boost relations for the quasilocal densities. A forth 
and final appendix presents a method for dealing with "corner" terms in gravitational actions 
(such terms are described below). 



II. PRELIMINARIES 
A. Foliations 

The boundary structure of M. leads to two classes of spacetime foliations. Our discussion 



of these foliations is close to one given by Hayward and Wong . 

Temporal foliations of Ai. The first type of break-up stems from a conventional ADM 
foliation of Ai into a family spacelike hypersurfaces. P A foliation of this class, referred to 
as a temporal foliation, is specified by a time function t : M. ^ I. The leaves of the foliation 
or slices are the level hypersurfaces of this time coordinate x° = t. Often, the possible time 
functions are restricted by the requirement that both t' and t" must be level hypersurfaces 
of coordinate time. The letter S is used both to denote a foliation of M. and to refer to 
a generic slice of this foliation, and the E slice specified by t = t^, (t* is some constant) is 
denoted E^^. If the manifolds t' and t" are level hypersurfaces of coordinate time, then it 
is convenient to set t' = Et/ and t" = T,t". The timelike, future-pointing, unit, hypersurface 
normal of a E foliation is denoted by u. 

Radial foliations of M.. The existence of the timelike boundary T suggests an alternative 
class of foliations of M.. Members of this alternative class are called radial foliations and 
rely on timelike hypersurfaces or sheets which have the topology of T (informally, sheets are 
radial leaves while slices are temporal leaves). One assumes that a radial coordinate x'^ = r 
parameterizes a nested family of such hypersurfaces which extend inward from T (see Figure 
(@)). This family of timelike sheets may converge on some degenerate sheet, and if this is 



In this paper we make a sharp distinction between the notion of quasilocal energy, which is sUce- 
dependent (observer-dependent) and unique only up to a reference point, and notions of quasilocal 
mass, which are boost-invariant and uniquely defined. However, both the quasilocal energy and 
the notions of quasilocal mass considered in this paper depend solely on gravitational Cauchy data. 
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the case, then the coordinate system breaks down at the degenerate sheet .0 With a notation 
similar to the one introduced above, we may represent T by 7^//, so the level hypersurface 
specified by r = r" is T (the inner radial sheet is The spacelike, outward-pointing, 

unit, T hypersurface normal is denoted by n (the unprimed letter n is reserved for a related 
but different vector field introduced below). 

Foliations of and T (see Figure (2)). It is of interest to examine how the S and T 
spacetime foliations mesh. If a temporal and a radial foliation of M. are simultaneously 
given, then the intersection i?f = f] is a two surface with the topology of B. Defining 
Bt = Bty, one finds that the family of Bt slices foliates T. This foliation of T and its generic 
leaf are both loosely referred to as B. The timelike, future-pointing, unit, hypersurface 
normal of this foliation is u. In general, the vector fields u and u do not coincide on 
T. Fixation of the time gives a family of sheets Br = Bt^^r which radially foliate the S 
hypersurface specified hj t = t^,. This foliation of S and its generic leaf are also represented 
by B. The spacelike, outward-pointing, unit normal of this foliation is denoted by n, and in 
general n and n do not coincide on E. 

Clamped foliations (see Figure (3)). Often in this paper we need to consider a particular 
subclass of S foliations, determined by the property that on T the timelike normal u is 
orthogonal to n. Such foliations are denoted by S with corresponding normal u. (So we 
have M ■ n = on T, where u is also the normal for the B foliation of T.) We described 
a S foliation as clamped. Note that it may not be possible for a temporal foliation to be 
clamped over all of T, since the u normals of t' and t" may not be orthogonal to n (assuming 
that t' and t" should be members of the family of S slices). We can also consider the locus 
of points which is the Eulerian history of B with respect to an (in general) undamped E 
foliation. This "boundary" , denoted by T is generated by the integral curves of u and may 
"crash into" or "emerge from" the actual boundary T. 

We maintain this barred and unbarred notation when it is necessary to deal with an 
undamped S foliation. However, in § III.B, which presents the derivation of the quasilocal 
densities, we make the clamping assumption which means that only clamped S foliations 
of spacetime Ai are considered (or every S foliation is a S foliation). When the clamping 
assumption is made, over-bars become redundant, and therefore in § III.B we drop all bars 
from the formalism. In this subsection of § III we assume that the u normals of t' and t" 
are orthogonal to the T normal n (in this section denoted simply by T and n). Though 
this is a limiting assumption, it in no way affects the generality of this paper, as we return 
to the fully general scenario in the following section. We demonstrate that the clamping 
assumption is a purely kinematical condition. 



^It should be emphasized that only a "local" radial foliation of an arbitrarily small spacetime 
region surrounding T is necessary for the analysis in this work. The full radial foliation of M is 
introduced only to have a closer analogy with temporal foliations. 
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B. Conventions and notation 



We adopt the following index notation. Lowercase Greek letters serve as Ai spacetime 
indices. Lowercase Latin indices from the latter half of the alphabet serve as S (and E) 
indices and as T (and T) indices. There is -hopefully- no confusion caused by this dual 
use of Latin indices. Lowercase Latin letters from the first half of the alphabet serve as B 
indices. Orthonormal (or when appropriate pseudo-orthonormal) labels and indices for each 
space are represented by the same letters with hats. For example, /i is a spacetime tetrad 
index and a is a. B dyad index. 

The spacetime metric is Qf^^, with associated (metric-compatible and torsion-free) covari- 
ant derivative operator V^, and e^*^ denotes a spacetime tetrad. The (pseudo) orthonormal 
symbol on spacetime is defined by eQjg-j = 1 = — e°^^^. Respectively, we have •jij and Vj 
{'jij and Vj), hij and Dj (Jiij and Dj), and aab and 6a denoting the metric and intrinsic 
covariant derivative operators on T (T), S (S), and B. We use ^^-^ {^f^'), Ef.^ (EfJ), and 
9a respectively, to represent a triad on T (T), a triad on S (S), and a dyad on B. Respec- 
tively, the permutation symbols on T (T), S (S), and B are defined by egj^ = — 1 = e^^^ 

(e6i2 = -1 = 6123 = 1 = e^'' (ei23 = 1 = and e^^ = 1 = (In Ref. § the 

convention for the T orthonormal symbol differs by a sign.) 



C. Spacetime decompositions 

The foliations just discussed lead to two decompositions of the spacetime metric. We 
choose to examine the metric in a coframe which has dt and dr as two of the coframe legs. To 
begin with, the temporal foliation S allows the matrix of metric components to be written 

as 

where the S indices run over (l,2,r). The N and are the ordinary ADM lapse and 
shift. Further, since each of the E slices is foliated independently by nested sheets with the 
topology of B, the matrix form of the E three-metric is given by 

Here, a and (3°" are the "lapse" and "shift" associated with the induced radial foliation of E. 
The super matrix formed by combining these expressions gives the so-called (1 -|- 2) -|- 1 form 
of the metric. The (1 -|- 2) indicates that three-space has been split into a radial direction 
plus a two-space, while the 1 indicates the time direction. 

Similarly, beginning with the full radial foliation T of spacetime one has. 
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where the T indices run over (t, 1, 2). The a and [3" are the gauge variables associated with 
this fohation. The submatrix associated with is 

where N and V"" are the lapse and shift associated with the induced B foliation of T . The 
super matrix of components for this splitting is the metric in 1 + (2 + 1) form. 

It is a straightforward exercise to express the "barred" variables in terms of the "un- 
barred" variables by simply equating the components of the (1 + 2) + 1 and 1 + (2 + 1) 
versions of the spacetime metric. First, define 

_V-n aV^ _ B-u NP' 

and the point-dependent boost factor 7 = (1 — f ^)~^/^ = (1 — -u^)"^/^. With this boost factor 
the set of transformation equations may be written as 

7V = ^ 

7 

a = 0:7 

V'' = V'' + V'P^ (2.6) 



(try) 
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The clamping assumption is tantamount to the f ^ (7 — 1) limit on the boundary T, 
in which case there is no longer a distinction between barred and unbarred variables. Note 
that in this case = V ■ n = aV^ = 0, which, as described in [0, implies that in the 
canonical form of the theory the S Hamiltonian can not drive field configurations across the 
boundary B. 



III. QUASILOCAL STRESS-ENERGY-MOMENTUM DENSITIES 

A. Action and variational principle 

Before turning to the derivation of the quasilocal densities, we must describe the action 
principle which is the cornerstone of our approach. Our starting point is the first-order 
Goldberg action [§] 



^ . r^Ae^Aa^, (3.1) 



where k = Svr (in units with G = c = 1) and F^.^^ = e^u ^ represent the spacetime 
connection one-forms which specify the Levi-Civita connection on M. with respect to the 
tetrad e^^. Also, the Sparling two-forms [y,|T^] are defined by 
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ap = --ep^f^r^^Ae^ (3.2) 

Therefore, as mentioned, the Goldberg action is not a Palatini action in which tetrad e^^ 
and connection r°" ^ ^ are varied independently. As it stands, the action ([Ol) possesses 
superfluous tetrad dependence. However, note that the Goldberg action is invariant under 
spacetime diffeomorphisms which preserve the boundary, since it is written purely in the 
language of differential forms. 



Our goal is to identify the Goldberg action ( p.l|) with the familiar "TrA'" action [0 used 
in metric gravity. The extrinsic curvature tensor associated with the S foliation is deflned 
by Kfj_y = —h^Vxu^j, (with the projection operator = + u^u^), while the extrinsic 
curvature tensor associated with the T foliation is deflned by B^i, = —7^ Va Ufj. (with the 
projection operator 7^ = (7^ — n^n^). The flrst step towards the desired identiflcation is to 
note that the action differs from the ordinary Hilbert action by a pure divergence 

r^Ae^Aa^ = ^/ ^e* - ^ f d(e^Aa^, (3.3) 
2k Jm 2k Jm 2k Jm ^ ' 

where e* is volume form on M.. Evidently, all of the action's tetrad dependence resides 
exclusively in boundary terms, 

- ^ / d (e^ A a^) = - / d^Xv/^V^ U^^ Va • (3.4) 

2k Jm ^ ' K JM ^ ' 

Now, if the time leg of the tetrad eg coincides with the future-pointing normals u on both 
t' and t", then the boundary terms associated with these manifolds are the desired TvK 
terms. Likewise, enforcing the condition that the third tetrad leg eg coincides with the T 
normal n on T ensures that the one obtains the desired TrG) term for the T boundary term. 
We assume that the variational set of tetrads obey these conditions. However, in general 
such tetrads are doubled- valued on the corners B" = t" ClT and B' = t' f] since u ■ n 
need not vanish on these two-surfaces. Therefore, to express the action (|3.1|) in the desired 
form, relax the second gauge condition on eg on a "small" (not connected) neighborhood 
of the corners such that the tetrad is single-valued. Next, take the limit that this small 
neighborhood "shrinks" to just the corners B' and B". Such a limit procedure is described 
in Appendix D, and it yields the following expression for the action: 

S^ = —f d^xy/^^+- f' d^xVhK f d^xy^Q-- r d^Xy/^(P, (3.5) 
2k Jm k Jt' K Jt k Jb' 

where <f) = tanh~^ v is the point-dependent boost parameter on B" and B' associated with 
the boost velocity v deflned in the last section. The corner terms were flrst given by Hay- 
ward and Wong for the metric action. Heuristically, they arise because, though the 



corners constitute a set of measure zero in the TrK integration over all of dA4, the trace 
of the extrinsic curvature is inflnite on these two-surfaces (as the normal of dAi changes 
discontinuously from m to n) . Note that the corner contributions to the action vanish if the 
initial and flnal slices are clamped to T. To obtain the variation of ( |3.5| ), one may straight- 
forwardly vary the action ( p.l|) and then apply the limiting procedure. This direct method 
is sketched in Appendix D. However, in the interest of brevity we borrow from results given 
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in Refs. Subject to the chosen "internal" gauge fixing, the action is a tetrad 

version of the metric action used in Ref. to define quasilocal stress-energy-momentum 
in general relativity. Hence, for the moment we may regard it as a metric action. Indeed, 
only the gauge- invariant quantities 7jj, h[p and h'l- are fixed on the boundary dj^ in the 
associated variational principle. Refs. P,p!^ have shown that the boundary contributions to 
the variation of are 



(58^) = r d^xp'^6hi,+ I d^xr^'57i,- - - r d\0(5v^, (3.6) 
V / dM Jt' Jr kJb' 



rt" r 1 rB" 

/ d=^xr^'57.,-- / 
It k Jb' 

where the gravitational momenta are given by 

^ 2k V 



(3.7) 



The variable p*-' becomes the standard ADM momenta in the canonical form of metric 
gravity, and it is conjugate to hij. Likewise, vf*-' is the ADM-type momenta conjugate to 
jij, but now canonical conjugacy is defined with respect to T. Note that equation ( p.6|) 
includes corner contributions to the variation which feature fixation of intrinsic geometry, 
in harmony with the fact that the induced metric is fixed on dAi. 

There is a complex- valued action functional, closely related to ( ^.1| ), which is based on 
the self-dual {+) (or anti- self- dual (— )) connection forms 

r(±)p- = 1 fYP- ^ ^eP^n^Yrf,) . (3.8) 



2 V 2 

This action is referred to as the complex Goldberg action and has the form (here we take the 
self-dual case) 



S 



^ . r%Ae7AaW^, (3.9) 
2k Jm 



where the complex Sparling two-forms are 

a(±^ = -ep^fAr^^^"'Ae^ (3.10) 



The complex action ( |3.9| ) differs from the previous one ( |3.1| ) by a purely imaginary boundary 
term. Indeed, setting 



S = S^-S\ (3.11) 



we find that 



S' = — d 
2k jm 



(3.12) 



With the gauge choices made above and the limiting procedure described in Appendix D, an 
appeal to Stokes' theorem yields 
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-S' = -^ / d'xVhe''^UspjE,^ + ^ _d'a:^r'Hspj^r\ (3.13) 

where frsj = C,fk '^j and Ursj = Dj Eg ^ are respectively the triad connection coeffi- 
cients on T and S. Notice that —S'^ contributes no corner terms to the action and that it 
serves as a subtraction term (a functional of the fixed boundary data) in the broadest 
sense (it depends on the boundary data of T, t', and t")-Q Because of the triad dependence of 
the subtraction term, we do not have the option of viewing the action (|3.9|) as solely a metric 
action. Furthermore, in order to fully remove the superfluous tetrad dependence associated 
with the action S, one would have to completely specify the triad on each boundary element 
of dA4 (though we do not chose to completely do so). 

Now consider the boundary-term contributions to the variation of the action ( p.9|) . Since 
the plan is to work with the Ashtekar variables in the canonical form of the theory, first 
express the boundary-term contributions ( p. 61 ) to the variation of the action in terms 
of the densitized triads on T, t', and t". (This is easily done with the identity ( [4.9|) given 
below and a similar identity for the T metric and triad.) Adding this result to the variation 
of ( |3.13|) , we find that 



{6S)gj^ = i r d^xA^, 6 [VhEr') + i /_ d'^xA'j 6 (v^^f - - T d\05v^, 

(3.14) 

where we have introduced the connections 

A' , = ^ (^%- - i ,) = i (-1 e''f u^gp, - i K'- ,) 

(3.15) 



X , = ^ (f + i ,) = 1(1 e''^ f,,, + i e 
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(With these conventions uj'^ gj = e*" pS j and f*" gj = e ps j-) The connection variable A^ j 
is (up to a factor of k) the S Sen connection, which becomes the Ashtekar connection in 
the canonical form of the theory. Likewise, the second connection j is the Sen connection 
associated with T. It is a complexified 5*0(2, 1) connection and enjoy properties completely 
analogous to the well-known ones enjoyed by the E Sen connection. In particular, in terms of 
the curvature of A" j one may compactly express the constraints associated the embedding 
of T in the Einstein space Ai. 0,|18[ Note that here these connections are not the canonical 
Ashtekar connections. We have not written down imaginary contributions to the corner 
terms which presumably arise from integration by parts on 6S^ terms. In fact, these vanish, 
and a calculation which demonstrates this is outlined in Appendix D. 



^To avoid confusion, it is crucial to note that in Refs. the notation represents an arbitrary 
subtraction term, while in this paper S'^ represents the specific term ( |3.13| ). 
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B. Quasilocal densities 



We now present all of the fundamental B tensors which serve as quasilocal densities 
describing the stress-energy-momentum content of the S gravitational fields contained within 
B. We express these densities in terms of the S Sen connection and triad. In the next section 
when studying the canonical form of the action principle, we consider the canonical versions 
of these expressions which are written in terms of the S Ashtekar variables. To begin 
with, we collect a set {s,ja, s"''} of quasilocal densities which is essentially the same as that 
described extensively in the original Ref. 0. This set is comprised of an energy surface 
density e, a tangential momentum surface density ja, and a spatial stress surface density 
s"^. We also find the need to introduce a new set {jh, Ja,t°'^} of quasilocal densities (also 
considered in [^]), which is comprised of a normal momentum surface density a tangential 
momentum surface density ja (which turns out to be the same as ja), and a temporal stress 
surface density t"'''. Both sets may be derived from the gravitational action ( |3.y| ) via the 
Hamilton- Jacobi method as described in the introduction. Therefore, we adopt the unifying 
point of view that any quasilocal stress- energy-momentum quantity is given by the rate of 
change of the classical action Sd corresponding to some variation in the fixed boundary data 
of DM = t'[j t" U T. However, we do not explicitly consider the classical action as in Ref. 0, 
since we prefer to "read off' the geometric expressions for above densities from the boundary 
contributions ( p.l4|) to the variation of S. 

In order to "read off' the various quasilocal densities from the boundary terms ( ^.14p , 
we make two assumptions in this subsection, (i) First, we assume that the E foliation of A4 
is clamped, so that u ■ n = 0. Again, this means that one may drop all overbars associated 
with three-boundary quantities from the formalism. Also, this sets = on the corners. 
The clamping assumption is made in this section only for convenience, and we return to the 
general slicing scenario in the next section, (ii) Second, we enforce partial gauge fixation 
of the triads on the boundary elements of Ai. Following Ref. we require that the T 
triad is time-gauge. This condition ensures that the T piece S'^\r of the subtraction term 
is functionally linear in the lapse N and shift V"'. As described in detail in Refs. [0,^ this 
linearity condition is crucial, because it ensures that the quasilocal energy density e and 
momentum density ja depend solely on the Cauchy data of S. Similarly, the triads on both 
t' and t" are required to be "radial- gauge." Essentially this just requires to coincide with 
n dX B. These restrictions on the T, t', and t" triads ensure that the purely imaginary piece 
of the corner contribution to the variation ( p. 14 ) vanishes (indeed we have already seen that 
this is a condition which follows from how the tetrad has been selected), and they ensure 
that the quasilocal densities to be defined behave appropriately under boosts. These points 
become clear below. The time-gauge and radial-gauge conditions are defined and discussed 
in Appendix A. Unlike the clamping assumption (i), these boundary gauge restrictions (ii) 
are absolutely necessary for our formalism. Once we have obtained both the geometric form 
and a physical interpretation of each quasilocal density, we turn in the next section to the 
issue of how these densities behave under boosts and also consider the canonical form of the 
action principle when the E slicing need not be clamped. 

Let us first examine the T contribution to the variation of the complex Goldberg action 
with the assumption of a clamped E slicing. Subject to the time-gauge condition, the T 
triad and cotriad can be expressed (at least locally) in terms oi a. B dyad 6a ^ and codyad 
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Oa 
y b, 



= i/N {d/dt - V d/dx"") 



Ndt 



C = 0\(dx'' + V^dt) . 



(3.16) 



The time-gauge condition has been indicated by replacing the triad label with _L. The 
associated time-gauge T connection coefficients are the following: 

rc±± = e^ [logiV] 



(3.17) 



ri2± = UN [aM 0[i '02]'- 1/2 e'^^ - V r^^, 

where for this set the "dot" represents partial time differentiation. Plugging these coefficients 
into S'^\rj-, one can verify that S'^\^ is functionally linear in the shift V"" and has no 
dependence. Next, applying the identities [rjfg = rf^ = diag{—l, 1, 1)] 



(3.18) 



to the T piece of the boundary variation (|3.14| ), we write the T contribution as 

N 



{6S) 



T 



J d^x^/a 



e6N-j, 6V'^--s'^' m 



ab 



Here the quasilocal density s"* is defined with respect to 

1 ss 



Sr 



-7 69^ 



T 



(3.19) 



(3.20) 



via s 



ab 



Sc 0^, and the expression {S6)ab is shorthand for 26aa5d"'b- Notice that 
{.5G)(ab) = ^cTafe, while {56)]^ab] is a pure gauge variation of the B dyad. Also, note that 
g[ab\ jg completely determined by the subtraction term 5S "\r (s("^) is determined by 5S^ and 
55*" contributions). Explicitly we have 



l_5S_ 



r 



13 



3a 



1 SS 



(3.21) 



r 



6S 



-1 



69^ 



T 



a"* - A' 



da n b 



We can rewrite these densities in terms of the S Sen connection. The appendix results ( P12|) 
express the time-gauge T Sen connection A^ j in terms of the radial-gauge S Sen connection 
A'^ j and other gauge variables. Insertion of the appendix results into the above expressions 
gives the following: 



Ja 



i# 



(3.22) 



= e^' A-ad a"- 6, " + {2i/K T^+^ - A,, 6, ' 



a 



ab 



Henceforth, we assume that e, ja, and s represent these expressions. Notice that e and 
ia are built exclusively from S Cauchy data {Ef^ , j)- Because of this fact, e and ja 
can be interpreted as canonical expressions depending on the Ashtekar variables. Because 
of the presence of 2ir*^"'") = % n-' + iTig^; the density does not depend solely on E 



Cauchy data. This term contains the spacetime acceleration = u" of as well as 

the T connection coefficient T|2x, which describes the rotation of the B dyad under parallel 
transport along the integral curves of u. Both of these terms depend on how the Cauchy 
data evolve in time. The real parts of the densities in the above set correspond exactly to 
the quasilocal densities first introduced in Ref. p|. Indeed, expressed in full detail. 



k 



K 



Vh 



rii (Taj P ■ 



■^i2a 



(3.23) 



ab 



1 
K 



ab 



be 



ab 



where kab is the extrinsic curvature of B as embedded in S (with k = a"-^ kab) and lab is the 
extrinsic curvature of B as embedded in T. 

We assume that each density has the same physical interpretation as given in Ref. 
and review these interpretations now. (For the following interpretations to be valid, one 
should consider the densities e, ja, and s"^ to be evaluated "on-shell", i. e. evaluated on 
some particular solution of the Einstein field equations.) From its definition ^/a e equals 
minus the time rate of change of the action S", where the time separation between the B slices 
of T is controlled by the lapse N onT (fixed as boundary data in the variational principle). 
Therefore, e is interpreted as an energy surface density for the system as measured by the 
Eulerian observers of S at B. The total quasilocal energy associated with the S gravitational 
fields is 
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E = d'^x y/ae , 

JB 



(3.24) 



the integral of the quasilocal energy density over the two-surface B. This energy is the value 
of the on-shell HamiltonianQ which corresponds to the choice = 1 and = on B , and 
it is a functional on the gravitational phase space associated with S. In a similar fashion, 
ja is interpreted as a tangential-momentum surface density, and the integral 



(3.25) 



is a total tangential-momentum quantity associated with the S gravitational fields and the 
B vector field 0". When 0" is a rotational Killing field whose orbits are contained in 5, the 
real part of rigorously represents the total quasilocal S angular momentum. p[ On-shell, 
the integral is minus the value of the Hamiltonian which corresponds to the choice = 0, 

= 0, and V^a^ = 0° on the boundary. In the special case when 0" is a genuine rotational 
Killing field, the form of ja makes it tempting to identify the imaginary part of with the 
"spin" of the B dyad. [|19| Finally, note that the real part of s"* represents the flux of the a 
component of momentum in the b direction. 

The original set of quasilocal densities have been obtained from a careful analysis of the 
T contribution to the variation (|3.14|) of the action. In a similar fashion we now analyze 
the t' and t" contributions to the variation.[] Often we drop the ' and " notations with the 
understanding that all expressions may refer to either the manifold t' or t". Remember that 
the t' and t" triads are radial-gauge. The radial gauge is indicated by replacing the triad 
label 3 by h. Therefore, with may split Ef.^ into a and (the gauge variables associated 
with the 1 + 2 split of hij) as well as the B codyad 6"" f With this assumption, we find 
identities like those in ( p.l8|) . Therefore, it easy to write the t' and t" contributions to the 
variation ( p. 141) as 



(55) 



d^^ 



jV 5a + ja5p''~-t 



ab 



mab 

The new quasilocal densities described at the beginning of this subsection are then 

1 6S 



(3.26) 



Ja 



±ab 



1 5S 



t" 



±0^a SS 



-iA\ea' 



(3.27) 



^Whether or not it is possible to find a truly satisfactory Hamiltonian for a spatially bounded 
slice S is a subtle issue in its own right. Following Ref. Q, this paper assumes that the correct 
Hamiltonian for a bounded region is the one which is "read off' from the canonical form of the 
gravitational action appropriate for a spatially bounded spacetime region. 

^The remainder of this section is based on Refs. |3,20|. 
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with the same expressions for the densities associated with the manifold t' . In full detail 
these are 



2 



UiaajP'^ - -^i2a (3.28) 



One could also write j\- = — where / is the trace of lab- Note that the definitions of ja 
and ja are identical, and hence ja carries the same physical interpretation as ja- Therefore, 
from now on we suppress the hat on ja- This equivalence of ja with ja results from the 
chosen gauge conditions. Also a result of these conditions is the fact that both jV and e are 
real. It turns out that the reality of j\- and e (or equivalently that the subtraction term 5''' 
has no a or dependence) is quite crucial, as it ensures that j\- and e behave well under 
boosts. Note that even if we had not enforced the radial-gauge condition on the t' and t" 
triads, then all of the densities listed immediately above would still by construction depend 
only on S Cauchy data. As shown in [Q], j\- is a normal momentum density, and the total 
normal momentum associated with the S fields is given by 



Jk = / d^x v^jV • (3.29) 

B 



This expression is minus the value of the on-shell Hamiltonian which corresponds to the 
choice A^ = 0, V^a'^ = 0, and = a l^'' = 1 on i? (heuristically, we may think of J\- as an 
on-shell value of the Hamiltonian which generates unit dilations of the system) . Finally, we 
refer to t*^* as the temporal stress density but lack a precise physical interpretation for this 
density. 



IV. BOOSTED DENSITIES AND THE CANONICAL ACTION 

A. Boost relations and invariants 

We now demonstrate that our collection of quasilocal densities behave under generalized 
boosts in a manner which is in accord with the equivalence principle. Fix a spacelike two- 
surface B in spacetime and also consider an arbitrary spacelike hypersurface E which has 
boundary (9E = B. The hypersurface normal of S is u. If we view the S slice as a member of 
a temporal foliation, then we may define the Eulerian history of B as T. By construction E 
is clamped to T. The observers at B who are instantaneously at rest in the S slice (Eulerian 
observers of S) determine the following set of quasilocal densities: 

e = 6°''^ Aafe 9c ^ 

j^, = -iA%ea' 
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Ja = iA^\ (4.1) 



t"^ = i (v4^ j Ef ^ o"^ - , a'^" Oa ' 



The primes appear on some labels in these formulae because in the triad formalism we set 
= u and E^i = n. | Note that here 2irj2X' = ajn^ + Itj^x^- Now consider a different 
hypersurface S which spans B (so like before dH = B). We may view S as a particular leaf 
of a temporal foliation which is not clamped to T, the Eulerian history of B with respect 
to E. Geometrically, the scenario now is identical to the bounded spacetime region M. that 
we have considered in the preliminary section. The observers at B who are at rest in the S 
hypersurface determine the set of quasilocal densities which are listed in ( p.22| ) and (p.27|) 



(simply the "unbarred" versions of the expressions above). We seek the transformation 
rules between the "barred" and "unbarred" densities, or, in other words, the behavior of the 
quasilocal expressions under switches of the hypersurface spanning B. With the appendix 
results ( |U14|) and ( |C15| ), it is quite a simple matter to establish that 



Ja=Ja--SaV (4.2) 

5«& = js''^-vje^ + - a"^ 7=^ u[v\ + - a"^ v 7^ n[v\ 

Kj Kj 

^ab _ ^ ^ab _ ^ ^ ^afe _ i ^ab ^3 ^ [^1 _ 1 ^ab ^ ^3 n ^ 

These are precisely the Eulerian-Eulerian boost relations found in Ref. 0. Remarkably, the 
particular form of the subtraction term ( p.l2|) , subject to the chosen gauge fixation, does 
not affect the boost relations of the "bare" densities. It must be stressed that if the gauge 
conditions (time gauge on T, and radial gauge on t' and t") had not been enforced when 



®At this point, the notation seems overly cluttered, but its use pays off in the appendix where 
the use of marked indices streamlines some derivations. In our notation spacetime quantities like 
r'^+) i2±' '^^^ barred. As the prime indicates, the T^^^ in (4.1) and the T^^^ in ( |3.22| ) are 
associated with two different spacetime tetrads. The labels i and 2 on the r*^^-* need not carry 
primes as the unprimed tetrad and primed tetrad share these two B legs, i. e. ej, = ej and eg/ = eg 
though in general e±i ^ e±. See the first two paragraphs of Appendix B for a fuller explanation of 
the notation. 
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defining the set of quasilocal densities, tlien tlie above boost relations would not have held. 
In particular, if e and jV are defined with subtraction-term contributions (which in this 
paper means they would no longer be real), then the first two boost relations are modified. 
For an application of the first boost relation to the Schwarzschild geometry see Ref. W^ - 
Following Refs. and 0, define trace-free parts of s"* and t""^ , 



rf = s'^'' - 1 /2 co''^ = - ik"^ - - A; a"'' + i " e^" 

K \ 2 



(4.3) 



Notice that the shear stress rj"'^ depends only on S Cauchy data. Of course, 9 = s"" a = 
1 / n {2 a^n^ — k + 2i rjg^) depends on how and the B dyad are extended into the future. 
Similarly, though t""^ depends only on S Cauchy data, its trace = t°- a = 1/k (2 6^ m'^ + / — 
2ici;j2i-) depends on how n and the B dyad are extended into the interior of B. We refer to 
as the shear temporal stress. Each of the densities e, jV, ja, v""^ and depend only on 
the extrinsic and intrinsic geometry of B and the normal at B^ One can easily see that 



(4.4) 

With the set {e^, jV, Ja, C°*} of quasilocal densities we can construct several invariants. 
For instance, notice that under boosts the density ja transforms like a gauge potential, since 
7^ 5aV = 5a4>- Therefore, the "field strength" or curvature Fah = 2 5\^a3h] of jh is an invariant. 
| 2T| Borrowing the results from P,|20[|, we write down the following list of simple invariants 



= - 

{m,f = VabV^'-CabC'' (4.5) 
{m2f = 2ie'''r]\C,, = 1/2 - (m,)^ . 

We make no claim that any of these invariants need be positive. One can also construct 
quartic-type invariants. In the metric formalism where all of the quasilocal densities are 
real, m^, (m^)^, and (mg)^ are linearly independent. For our theory, the real invariant 
along with the real and imaginary pieces of (m^)^ comprise a linearly independent set. It is 
natural to add to this list the scalar curvature TZ of B. One invariant of interest which we 
can build is 



(^^^ C^A<xK = fc"'' kab -k^ - lab + l^ + TZ = -^^2 + Reimsf - 2k-^ u] , (4.6) 



thank J. D. Brown for making this point. 
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where C^xa^ is the Weyl tensor of g^^ (here we are assuming vacuum), the two- metric 
^^lu = 9fiu — n^rii, + here serves as the projection operator into the B shces, and Re 



stands for "real part." This equation is a geometric identity associated with the embedding 



of the two-surface B in spacetime Ai. [p| -p3| As described in Ref. ||22|, the Bondi mass 
arises as the appropriate null asymptotic limit of the quantity 



M 



A 

16^ 



(4.7) 



The factor A is the area of B, and it has been inserted in order that the above expression 
has units of energy. Moreover, since the invariant = (8/K^)/ip is essentially the product 
of the convergences p and —fx of the null normals {u ± n)/^/2 to B, another invariant we 
can built is the Hawking mass [p!3| , p2 



A 



H 



IGvr Jb 



(4.8) 



The Hawking mass also yields the Bondi mass in a suitable limit. |22 



B. Canonical action 



Our goal in this subsection is to consider the variational principle associated with the 
canonical form of the action S ( p.ll| ). In order to express S in canonical form, we first 
consider the (3+1) form of S. Begin by expressing in (3+1) form. This can be done 
straightforwardly with standard methods. However, a faster way to get the (3+1) form is to 
borrow the results from Ref. P, p!^ , p!8[] . In those references the action 5*^, viewed as a metric 
action, has been expressed in canonical form.^° Therefore, we may simply cast this result 
into the language of triads. This is achieved by assuming that the S metric is a secondary 
quantity derived from Ef ^ and by using the identity 



dhij/dEf ^ = (E) ^{hij E^ k — hkj E^ j — hik E^ 



In anycase, the result is 
1 



d*x 



M 



K 



K"^ jEf^ - NH-V^Hj 



+ J^d^x 



(4.9) 



(4.10) 



where we have the following: 

n 



1 

'2k 



^'^For the metric action the (3+1) form is obtained from the canonical form by simply assuming 
that p^^ has the form given in (|]^) and that 2NKij = —hij + DiVj + Dj Vi. 
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(4.11) 



Hi 



Here R is the Ricci scalar of E, and , and {j^)b stand for the real parts of the 

densities in ( p.23| ) and ( |3.28| ) (the notation is redundant for e and as these are purely 
real). Also, at this stage, the hybrid extrinsic curvature j is merely a short-hand notation 

for a complicated function of A^, , Ef\ and E^K Finally, note that and are not 
constraints. 

The next step is to calculate the (3+1) form of the subtraction term, which can be 
written as 



-S" = -^ f d'^xd, {e^'P cusp, EfA + ^ [ d 
Tedious but straightforward manipulations yield 



V / ' spr ■ 



(4.12) 



2k Jm 



2k Jm \ J 2k Jt 



X \J ^ 6 '^s'pf 1 



(4.13) 



where = h E^^ . For the middle integral on the right-hand side, we now use Stokes' 
theorem for each S slice and enforce the radial-gauge condition at the boundary B of each 
S. Also, we expand the integrand in the final integral subject to the assumption that the 
T triad is time-gauge with respect to u. The result of these calculations is 



S" = -— d^xe^'Pusp^E,^ 
2k Jm 



2k Jr 2k Jt 

Next, "barring" the last formula in (|3.17|) , one finds 



atd ea 'e^' = N f^^^, + 1/2 1'^^ 6aV, + V f 



12b 



(4.14) 



(4.15) 



Insertion of this formula into ( |4.14] ) gives the desired (3+1) form of the subtraction term, 

-S° = -— I d^x e^'P uJsr,, Ef^ ~- I d^x v^c^tofc , (4-16) 
2k Jm k Jt 

where we have also used fact = ^acb- 

We now turn to the canonical form of the action principle. We shall avoid the issue of 
the reality conditions by working first with canonical form of the real action . Therefore, 
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upon adding the pure imaginary boundary term ( [4.16| ) to , we merely introduce a complex 
chart on the real phase space. [Q The canonical form of the action is the following: 



d^x 



M 



P' jEf^ - NH-V^ Hj - 1 /2 r 



In general j ^ l/nK^ j. Indeed, setting Pij = P^ j Ef^i, one has that 

P ^ -K 

Hi 

where ~ stands for modulo the rotation constraint 



jfs = 2P\^- -5'^^ Ei^ 



(4.17) 



(4.18) 



(4.19) 



(Note that the Lagrange parameter (pfg associated with J*"** is not related to the boost 
parameter and that 5"''^ is the Kronecker symbol.) Furthermore, in ( 4.171 ) T^i 'Hj, 7i^, and 
Til have the same forms as given in ( |4.11| ) but now are built with P^ j rather than j. 
In particular, in the canonical action 



Hi 



— ■yk + vy aij P*-' 



rii aaj P'^ 5a(p 



(4.20) 



(with the radial- gauge condition at the boundary S of S, one can write k = —uj"'\-a). 
As is well-known, the anticommuting Lagrange multiplier (pfs associated with the rotation 
constraint can be geometrically interpreted as the time component of the connection forms, 

m 



-r,,t = - (r,, , d/dt) 



(4.21) 



Enforcement of the radial-gauge condition at the boundary B of each E slice places a bound- 
ary condition on (pfg- This boundary condition is the canonical version of setting the con- 
nection coefficient Tah/iU^ = 0, where we are working in the RT-gauge described in the 
appendix. The coefficient Ta\-fj.u^ describes the rotation of n as it is parallel transported 
along the integral curves of u. To see what the required boundary condition on (pf^g is, first 
recall that in the triad formalism the vector constraint is not the generator of diffeomor- 
phisms, rather Ti.f'^ = Hj — 1/2 J^'^ cursj so the vector constraint generates rotation of the 
triad. |^ Therefore, the boundary condition 



- 1/N {^(j),^ + u;,^,V') 



r 







(4.22) 



ensures consistency between the selection of the radial-gauge condition for the S triad at B 
and the evolution of the triad as obtained from the variation of the canonical action. 
We now add the boundary term ( 14.16| ) to the canonical action ( [4. 17] ) and get 
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S = [ d^x 
Jm 



+ J d^x 



K 



^ -NH-V^Ub 



(4.23) 



where in anticipation of dealing with the Ashtekar versions of the canonical constraints, we 
have written N = h'^/^ N and H = h^^^H. Here H = , while Hb = Hi + i/K^LUi^h- 

Furthermore, for the rest of this section j is the canonical Ashtekar connection 

1 



A 



K 



As usual, one may replace the rotation constraint with the Gauss constraint, 



K 



2 ^rst ) 



j'st 



(4.24) 



(4.25) 



where "^Vj is the derivative operator associated with the Ashtekar connection. Moreover, 
using the Ashtekar curvature. 



tk 



(4.26) 



2D[jUj'^ fc] + e'"''* ^sj ^ik 



sj ^ ik ^ i ^[j k] ; 

one can build the standard S constraints: 

C = 1 e-* Er ' % ' =n-iD, [Er ' f) 



(4.27) 



With this machinery, one may rearrange terms in the expression (|4.23|) to find 



S = J^d'^x^iA'jEf^ - NC-V^Cj-ip''Jfj + y_d^x 



--v^- NC-V''Cb 



(4.28) 



The Lagrange multiplier associated with the Gauss constraint here has the explicit form 

(f^ = -1/2 e^^P (f)sp - i 5'''Es ^ Dj N + ikP^ j . (4.29) 

Furthermore, now we have 

C = Vhn-i-f{dr)ief'PP^jEs^Ep' = aa {-f e - v-f j^) 



Ca ='Ha = —\f(y 



ja Sa(, 

K 



where a = h^^^ a (Again, C and Cb are not constraints, i. e. they do not vanish "on- 
shell.") At this point the densities e, and ja have the same forms as in ( p.22| ) and (p.27| ) 
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but are constructed with the canonical Ashtekar connection. Therefore, off the constraint 
surface in phase space defined by the Gauss constraint, the energy density e is no longer 
manifestly real. Notice that C has been defined with a higher density weight, because it 
is paired with the boundary "Lagrange multiplier" N, which we have taken as a density 

of weight minus one. We also remark that the kinematical torsion which is present in the 
Ashtekar connection modifies the boost relations. Therefore, for instance, it is not true that 
e = •y e — vy j\- in the canonical picture. 

Before considering the variation of the action ( 4.28| ), we find it convenient to rewrite the 
Lagrange parameter in the following way. Take 

(f)fs = -Tfst = -NTfs± - i^fsj , (4.30) 

and also write 

- i ^ Dj N = -iNa^ = -i Nr x± , (4.31) 

where af = Ef [log N] are the triad components of the spacetime acceleration of u. With 
these relations one can set 

^ /-sp -kA^j , (4.32) 

which is, of course, essentially the well-known result that <y9^ = e'^'^^T^'^^ gpt- We shall need 
the expression for (p^ when the radial gauge condition is enforced, 

if^ = 2A^r(+) -f^A^j ■ (4.33) 

Using = + V'^ al, one can put this result in the handy form 

ip^ = 2iVr(+) - Ef ^ + m (> + je V"") . (4.34) 

Direct calculation yields the following for the variation of the canonical action ( |4.28| ): 

5S = (terms which give the constraints and equations of motion) + i j d^x AC j SEf 
-J.d'^x C6N + Cb 6V^ - v^iV/2 (7 s"'' - vy t"^ + 1/k a"^ + A a"^) {56) 

+- f d^x ^/^\KNhj^-vye)+SaV^- 1/2 a^^ &ab\ S(f) d'^x(j)6^/^, (4.35) 

K JT ^ J K JB' 

where here N is h^^"^ N, and it would perhaps be better to express {S9)ab as a variation in 
terms of the densitized dyad (as is certainly possible). Also above, 

NA = iA^ J Ef ^ + jV + je + i/t^ (v^ - 2N T^+^ . (4.36) 



Let us verify that our result for the variation of the canonical action agrees with the variation 
( p.l4| ) of the non-canonical action. With the interpretation ([4.34|) A vanishes. Therefore, 
enforcing the Gauss constraint and using the boost relations ([4.2|), one finds that 
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{6S) 



T 



r 



d^x \e5N - jb 5V'' - N/2 r'' {5e)ah 



+- / d^'x^UN j^> + 6aV^- 1/2 a''''aab)H-- ^ d^x<j)6V^, (4.37) 



where now one must again consider the quasilocal densities to be expressed in terms of the 
Sen connection. The density j^i = — l/Kcr"^/ab, and in the non-canonical picture 



2NI 



ab 



-O-ab + SaVb + 5bVa ■ 



(4.38) 



so the middle integral on the right-hand side vanishes in this case. This means that (p is 
not held fixed in the canonical variational principle, as the equations of motion ensure that 
this term vanishes for arbitrary variations 6(f) about a classical solution. Using the barred 
versions of ( |3.21|) , we then have 



T 



K JB 



B" 



d'^X (j) , 



(4.39) 



in agreement with the variation ( |3.14| ) of the non-canonical action. 



V. DISCUSSION 

We conclude with (i) a description of some of the new developments in the theory of 
quasilocal stress-energy- momentum which will appear in Ref. [^. We also briefiy comment 
on several technical matters. These are (ii) the issue of additional subtraction-term contri- 
butions to the action, (iii) the relationship of our formalism with the Sparling two-forms, and 
(iv) a problem encountered in the attempt to extend the Brown- York notion of gravitational 
charge to the Ashtekar- variable version of the theory. 

(i) Since much of the analysis in this paper is based on Ref. it is appropriate to 
describe a few results which will be found in this upcoming work. First, Ref. deals 
exclusively with the metric-variable version of quasilocal stress-energy-momentum, though 
this is not a distinction between Ref. and the present paper that we wish to highlight in 
this paragraph. Regardless of the choice of gravitational variables, the results to be found 
in Ref. are more general than those presented here in the following sense. In this paper 
the set {e, jh, ja, s'^*, t"''} of quasilocal densities is associated with a family of Eulerian (or 
surface-forming) observers at B. We have derived the transformations rules between two 
different sets of quasilocal densities, but each of the sets is associated with its own family of 
Eulerian observers. Ref. 0] also considers a set of densities {e,j\-,ja, (but built with 

metric variables). However, in Ref. the densities need not be associated with a family of 
Eulerian observers. That is, they may describe the stress-energy-momentum content of the 
gravitational field which is associated with a family of Lagrangian (or non-surface forming) 
observers, such as those determined by the timelike Killing field of the Kerr geometry. The 
transformation rules between the Lagrangian set of densities and a set associated with an 
arbitrary family of Eulerian observers will be given. Hence, the boost relations which will 
appear in Ref. are more general than those appearing here (the Eulerian-Eulerian boost 
relations arise as a special case) . It may possible to derive these more general boost relations 
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in the Ashtekar-variable framework as well, however, such a derivation would be considerably 
more difficult from a technical standpoint. 

(ii) As mentioned in the introduction, we have chosen to formally treat the imaginary 
boundary term —S^ as a true subtraction term a la Brown and York. However, we now 
argue that in some contexts it is necessary to consider the freedom to append to the action 
S = — an additional subtraction term —S'^. In the interest of economy we restrict 
our argument to matters concerning the quasilocal energy surface density e, though perhaps 
much of this discussion also pertains to the other quasilocal densities. Consider ffist the 
Brown- York expression 



-(k-k"). (5.1) 



K 

In the metric formalism, as in this paper, k represents the trace of the extrinsic curvature 
of B as embedded in E and comes from a base action in the derivation. The k^ term 
represents the trace of the extrinsic curvature of a two-surface which has the same metric as 
B, but which is uniquely embedded in a three-dimensional manifold possibly different than 
S. In the Brown- York formalism it arises from a real subtraction-term contribution to the 
action. When possible, this real subtraction term is typically chosen such that the different 
three-space is R^, and hence the k'^ term references the energy against fiat-space. For a 
given asymptotically- fiat spacetime, the presence of the appropriate k^ term is crucial if the 
quasilocal energy, 

E = - f d^x^(k-k") , (5.2) 

K JB ^ ^ 



is to agree with the ADM notion of energy in the suitable limit. |pi)|,H;|25 

Though we have added the imaginary boundary term —5*^ to in this paper, the 
resulting quasilocal energy, 

E= f d^xy/^e^^AabOc^^- I d^xy/^k, (5.3) 

JB K JB 

is really only the "unreferenced" energy. (As we have seen, the particular subtraction term 
used in this work makes no contribution to e and thus E.) If we wish to put the Ashtekar- 
variable expression for the quasilocal energy into full accord with the ADM notion of energy, 
then we should allow for the freedom to append to the action yet another subtraction term 
—S'^. Use of the more general action S' = S'^ — S"^ = — — S'^ in our analysis would 
yield 

E = j^d^x^ e^'^ e, ' [Aab - A?,) (5.4) 

for the Ashtekar-variable quasilocal energy. The new reference-point contribution e"*^ 6^ ^ 
stems from —S'^. At this point —S"^ is an arbitrary functional of T data. With this new free- 
dom, the quasilocal energy can be defined to agree with the ADM energy for asymptotically- 
fiat spacetimes (in the suitable limit that B becomes the two-sphere at infinity). A fuller 
discussion of this issue will appear elsewhere. 
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(iii) The real and complex Sparling two-forms obey the Sparling relation 

da^ = daW/i = r^ + G'>e;, (5.5) 

where e% is a basis for three-forms, is the Einstein tensor, and are the Sparling three- 
forms. The explicit form for (which is real) is not needed here but may be found in, for 
example, Ref. The Sparling relation suggests that (when pulled-backed to a three- 
dimensional slice S of spacetime) may be interpreted as a tetrad-dependent local energy- 
momentum density for the gravitational field. ||r5|,|9[| The corresponding frame-dependent 
potential can be taken either as or a'-^-* ^. If we fix a two-surface B and its spanning 
three-slice S in spacetime, then the boundary structure of our selection provides a natural 
(almost-unique) frame at B. Namely, the radial time-gauge tetrad of Appendix A, which 
has the E hypersurface normal u as its time leg and ra, the normal of B in S, as its third 
space leg. With this frame choice, the pullbacks s*{a^^^ ^) to B (s is the inclusion mapping 
s : B ^ M) are the following: 

— KB \fa 

njhVo'd^x (5.6) 
kO^^ jb^/ad^x . 

Here these are expressed in terms of the S Sen connection and triad. Further, the pullbacks 
s*(cr^) of the real Sparling two- forms are just the real parts of the above expressions. But 
only the last expression is complex, so s*{a±) = s*{a^~^^ ±) and s*{a^) = s*{a^~^^ One has 
s*{o'a) = K,9a^{j^)b d'^x. Scc Rcf. for more details. 

(iv) The Brown- York notion of gravitational charge is based on the T momentum con- 
straint, 

- 2 A (r, - (7r^)%) = -l/^^^^n'G,x , (5.7) 

where vf*-' is given in ( p.7|) and, in the metric formalism, [n'^Y^ depends only on (and so it is 
annihilated by T>i). Now we work on-shell and in vacuum, so this expression vanishes. Brown 
and York define a "stress tensor" f*-' = 2/ a/— 7 (vf*-' — (vf'^)*-'). Assume that T possesses a 
Killing field , and so T>if^^ (j = 0. Therefore, since —Uif^^ = eu^ + jbO'^\ one has the 
following conserved charge: 

Q^{B) = I d^x^ [eu^ + jb o'') . (5.8) 

When attempting to introduce such a notion of charge into our formalism, we run into some 
difficulty since the subtraction term S"^ may be triad-dependent. (Here S'^ may or may 
not represent the particular subtraction term considered in this work.) The natural way 
around this difficulty is the following. First define 

{Wy, = 5S'"/5{^i,^) . (5.9) 

In our situation (n'')^^ = {Tl^Y j is not necessarily annihilated by Vj, though Vj (n^)(jfc) = 
0. Therefore, set {tt^Y^ = ^—7/2 (^{H'^)^ j.'y^^ — (H'^)*^*-'^^ and use it in the above construc- 
tion. The charge may now be complex, but, subject to the assumptions made above, it 
is conserved. 
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APPENDIX: KINEMATICAL FRAMEWORK 



Appendices A, B, and C outline a kinematical framework for examining how the intrinsic 
and extrinsic geometry of spacetime as foliated by a family T hypersurfaces is related to 
the intrinsic and extrinsic geometry of spacetime as foliated by a family of S hypersurfaces. 
With this framework one can express objects such as the T extrinsic curvature Bjj or 
the time-gauge T Sen connection j in terms of the intrinsic and extrinsic geometry of 
spacetime as foliated by S hypersurfaces. Such a "splitting" of the T Sen connection is 
needed in order to derive a similar splitting of the S radial-gauge Sen connection A^ j in 
terms of the geometry of the S foliation. The splitting of A^ j is used to obtain the boost 
laws ( [4. 21) for the quasilocal densities. The kinematical framework consists of (i) two distinct 
spacetime tetrads (one adapted to the S foliation and one adapted to the T foliation), (ii) the 
transformation equations between these tetrads, and (iii) the inhomogeneous transformation 
law between the sets of associated connection coefficients. The relevant spacetime tetrads 
are constructed in Appendix A, and their associated connection coefficients are tabulated in 
Appendix B. Appendix C outlines the splitting procedure by applying it to Qij, the simplest 
example. We then quote the splitting results for the T time-gauge connection coefficients 
ffsj, j, and A"^ j. The final Appendix D applies some of this formalism to explain the 
origin of the corner terms in the action 



APPENDIX A: ADAPTED TETRADS 



The boundary structure of Ai suggests two natural classes of spacetime tetrads. The 
first class is a subclass of time-gauge tetrads determined by the boundary structure of E. 
The second class is a subclass of "radial-gauge" tetrads determined by the B foliation of T. 
These tetrads need only be defined on some small spacetime neighborhood surrounding a 
portion of T. We do not address the issue of whether or not either of these tetrads can be 
extended globally over all of Ai. 



1. Radial time-gauge tetrads 

Enforcement of the time gauge condition locks the time leg of the tetrad to the S foliation 
normal u. This condition is indicated by replacing the tetrad time label with ± so that 
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e_L = u. Because each S shce has a boundary _B, a natural subclass of all time-gauge tetrads 
exists which is determined by an auxiliary condition on T . This further requirement is that 
on the three-boundary T one of the space legs of the tetrad, chosen to be e\- = eg, coincides 
with n. One should note that this correspondence is not made between eg and n in general. 
Such a choice of tetrad is said to obey the radial time-gauge or RT-gauge (see Figure (H)). 
RT-gauge indices and labels take the values (±, 1, 2, h). Now the usual assumption is that the 
vector field d/dt points everywhere tangent to the hypersheets of constant r. Equivalently, 
(dr, d/ dt) = or dr /dt = 0, and the r coordinate is Lie transported along the integral curves 
of the time vector field. This assumption results in almost no loss of physical generality. 
It does demand that the integral curves of the time vector field may not emerge from or 
flow into the three-boundary T. However, since the spacetime-filling extension of the three- 
boundary T in terms of hypersheets of constant r is completely arbitrary, on the interior of 
S these integral curves can be chosen to flow in any direction (as long as the tangent field 
d/dt lies at each point within the future light-cone). Subject to the requirement dr /dt = 0, 
one can write the most general radial vector field mapped to unity by dr as d/dr = an-\- f5, 
which is similar to the familiar formula d/dt = N u-\-V. As seen earlier, the variables a and 
are respectively the kinematical "lapse" and "shift" associated with the induced radial 
foliation of the S slices. Therefore, we have the following explicit formulae for the RT-gauge 
tetrad and cotetrad: 

^± = l^[^,-y'E,-V^E,^ e^ = Ndt 

ea = Ea = Oa c^ = 6^ + V^dt + /3^dr (Al) 

e^ = E^ = -(-^-l3^ea] = adr + V^dt. 

a \dr J 



2. Time radial-gauge tetrads 

The radial-gauge condition requires that one of the space legs of the tetrad, taken to be 
Ch' = eg , coincides with the T normal n. A natural further requirement can be placed on 
radial-gauge tetrads. Namely, the time leg e±r = eg can be tied to the B timelike normal u, 
so the indices and labels associated with this class of tetrads run over (±', 1, 2, h'). Such a 
tetrad is referred as time radial-gauge or TR-gauge (see Figure (^). Now the radial vector 
field is written as d/dr = an-\-(3, though it still points tangent to the S slices. The variables 
a and are associated with the T foliation of Ai. On T one can express the time vector 
field as d/dt = Nu + V, where N and V"' are the gauge variables associated with the B 
foliation of T. The RT-gauge tetrad and cotetrad is 

e±' = L' = ^ - = + /9^'dr 

ea = ^a = 0a 6^ = 6^ + V^dt + ^^dr (A2) 
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We explain the chosen notation further in the next two paragraphs. 



APPENDIX B: ASSOCIATED CONNECTION COEFFICIENTS 

For the special tetrads considered above, certain of the corresponding connection coeffi- 
cients have notable geometric meanings. This subsection is a glossary of various connection 
coefficients and their geometric interpretations. RT-gauge connection coefficients are rep- 
resented as r^CT^, while TR-gauge connection coefficients are represented by 5-/^/. Note 
that inspection of the indices allows one to discern which set of connection coefficients is 
being dealt with. 

Perhaps a few more comments on the notational scheme will be clarifying for the reader. 
The RT-gauge tetrad and TR-gauge tetrad e^/ are both tetrads on the same spacetime 
Ai. As has been evident, "e" is used to denote both tetrads, and it is the type of label 
(primed or unprimed) carried by "e" which makes the notational distinction between the 
two tetrads. Clearly then, though the 5-^ and 5-/^/ are different sets of connection 
coefficients, they specify the same spacetime connection (that of Levi Civita); and so we 
use r for both, again letting the labels make the notational distinction between the two 
sets. Notice that neither e-*" nor need be surface forming. For instance, A de-*" 7^ 
in general. Technically put, in general and will not satisfy the Frobenius condition. 
Therefore, e-*" is not necessarily of the form Ndi, where t is a coordinate which specifies 
some S slices. However, this is unimportant for our calculations. What is important, is 
that the clamped E slices have a future-pointing normal u which agrees with e±> at the 
physical boundary T. In fact, the E slices are determined by this condition, but clearly their 
extensions off the three-boundary are highly non-unique. Ultimately, we are interested in a 
S 3-slice and S 3-slice which span the same B 2-slice of the timelike boundary T. We want 
to compare the Cauchy data of the two slices at the mutually bounding two-surface B. Our 
formalism handles this issue but does not assume that the S normal u agrees with e±i off 
the boundary. 

In the lists of this appendix, since the geometry of M. is torsion-free (i. e. the torsion 
two- form of Cartan vanishes [1^), all of the extrinsic curvature tensors are symmetric. Note 



that for the extrinsic curvature tensors defined below, we adopt a different convention for 
the staggering of indices than the convention used in Refs. 0,0]. However, since all of these 
tensors are symmetric, all of our results match those found in Refs. 

1. RT-gauge connection coefficients 

The RT-gauge connection coefficients are tailored to B as embedded in S. We have the 
following correspondences: 
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(Bl) 

Note that the formulas for ^ and a*" are general time-gauge expressions. Also, are the 
tetrad components of the spacetime "acceleration" of n, while the S "acceleration" of n has 
components b"". For fe*" the f is a T index and can take the values (_L, a), while for s and 
a^' the f and s are S indices taking the values (1, 2, h). 



2. TR-gauge connection coefficients 



The TR-gauge connection coefficients are tailored to B as embedded in T. We have the 
following correspondences: 

(B2) 

t c = — e" u eg V,, e_L' = — _l'c = — t"" ±'c 



Like before, the formulas for 0'' § and b'^ are general radial-gauge expressions. For 6*" § and 
6'' in this list the f and s are T indices taking the values (_L',a). The (f can take the 
values (a, h')) are the tetrad components of the spacetime acceleration of u, while the T 
acceleration of u has components d'^ = ±'±'. 



APPENDIX C: SPLITTING PROCEDURE 



1. Transformation equations 



The set ( p.6|) of transformations for the metric variables can be used to express the 
transformations between the RT-gauge tetrad ( [AID and the TR-gauge tetrad ( |A2|) . For 
example. 



N \ dt dx^ 
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N \dt dx^ 

7 ,A 
N \dt dx^ 

7 ex + f 7 ei- . 



The complete set of transformations is 



e_L' 



7 e_L + f 7 Ch 
f 7 e_L + 7 Ch 



1^ 

dx^ 



or or ox" , 



■,-L' 



e = 7 e — f 7 e 
= — f 7 e-*- + 7 e'" 



(CI) 



(C2) 



Notice that the B legs of both the tetrads are the same, which is why the notation can be 
compressed so that TR-gauge tetrad indices like p' run over (±', a, h'). 

The inhomogeneous transformation rule describing the behavior of the spacetime con- 
nection coefficients under the above tetrad transformation is the following: 



(C3) 



This law provides the bridge between the TR-gauge connection coefficients (|B2|) and the 
RT-gauge connection coefficients ([Bll). 



2. Geometric link between T and S 

As an example, we apply the developed formalism and derive the splitting result for 
the three-boundary extrinsic curvature Bjj. This result has been obtained via ordinary 
tensor methods with projection operators in Ref. However, the ordinary projection- 
operator method is not sufficient for calculating the analogous split of the T Sen connection. 
We provide the splitting calculation for G)jj here as a simple demonstration of how such 
calculations are performed. Beginning with the first expression of ([B2D , one uses the rule 
( PSD in tandem with the set (|C^) to find 

I = -e% e, (v^ x^ + 7 ka) -e' ^ e, [v^] -e\e, [7] • (C4) 

(Note that in this equation f and s are T triad indices which take the values (_L', a).) A bit 
of work and the relations ( Pl| ) yield the set of B'' s triad components, 

G)"*" X.I = —7 + "^7 -^'^ h — 7^ ej_ [f ] — f 7"^ ey- \v\ 
e^'a = ir^-7'eaM (C5) 
e'c = 7^S + ^7^''c- 

With the set (|C5|), construction of the sought-for splitting of Bjj is not difficult. For conve- 
nience work in spacetime coordinates. The boundary three-metric may be written as 

li,v = CTf^u -U^Uy , (C6) 
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where the two-metric a^,y = g^y — n^riy + u„ here serves as the projection operator into 
the B shces. Wiring the above form of 7^, the identity operator on T, on each of the free 



indices of 6^,^, one obtains 



Q^,u = Uy ex'±' - 2 (T^) 0x' A + crjcr;: 0a« , (C7) 

where an appeal to the symmetry of 9^,^ has been made. Plugging = 7 + f 7 and 
the results from (|C5|) into (|C7| ), one arrives at the following split of the three-boundary 
extrinsic curvature: P^,|5[| 

Q^u = lkf,y + f 7 Kij cr^ al 

} (C8) 
X I7 — f 7 + 7^ u[f ] + f 7^ r;,[f ] | 

+2 {7 + t;7 } {n^ i^ij - 7^ A^^} 

Enforcement of the clamping condition f ^ recovers equation (A. 16) of Ref. 0, 

0a/3 = + Up + 2 (T^) n-' K.;^ . (C9) 

The set of T time-gauge connection coefficients is |f"_L/_L/ , _L/e , r'^c±' , ^"gs}' where 
the first two have been considered in the set The splittings of these expressions are 



< «P = 7^ < Op + {v-if hp 
l^u = 7 Krp crj^ crP + v-f kf,^ 

fdc±' = 7 TacX + VI (^ach (CIO) 
Tacb ^acb ■ 

Using this set and (|C8| ), one finds the following split of the time-gauge T Sen connection 
(pulled back to B) in terms of the radial-gauge S Sen connection (pulled back to B): 



(Cll) 



<^ A = -a;^ (A^ + (i77«^)VAt; 

a'^A' A = {v^A' a + i7e''^cA) • 

Now consider 2r*^+^ = ^i2± " i'^-'^, and the following clamped results (f — > limit) for 
the splitting: 



^^, = -«,(2/«:)rWi2^-aM' , 



^A AO, 1 „\^ac A 

c^llA a = -io-„e Acx 



^ A 

(C12) 



To find the splitting of the needed pieces of the radial-gauge S Sen connection in terms 
of the radial-gauge S Sen connection and other gauge variables, first find the split of A 
(pulled back to B) in terms of the T foliation variables, 



r 

j 
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(C13) 



Combination of this result with (|C11|) gives 

(C14) 



The boost relations (|4.2| ) for e, and ja can be derived with these expressions. To derive 
the boost results for s"* and t"'', one must use these expressions and also the results 

(C15) 

kA^' ^n^ = 7Ky4^ ^n'' - v'y2T'^'^^ + 17^^ M 

Note that on the left-hand side the selfdual coefficients are TR-gauge, while those on the 
right-hand side are RT-gauge. 



APPENDIX D: CORNER TERMS IN THE GRAVITATIONAL ACTION 

This appendix presents a simple tetrad method for analyzing "sharp-corner" terms in 
the gravitational action principle. We show how the corner terms in the action (|3.5|) arise. 
(Using a different method, Hayward has made a systematic study of such corner terms. 
As mentioned, the Goldberg action differs from the Hilbert action by the pure divergence 

d(e^Aa^)=-/ d^Xv^V^fe^'^e^^r^,) . (Dl) 

To ensure that, upon the use of Stokes' theorem, this divergence gives the desired ^^TrK" 
and "Tr0" terms on the boundary elements, tie 65 to the u = e\_ hypersurface normals 
on t' and t" and tie eg to the normal n = e^i on T. However, if these gauge conditions 
are enforced simultaneously, then in general the tetrad is doubled-valued on the corners 
B' = Bf^r" and B" = Bf^r"- Therefore, in order to both retain the desired ^^TrK" and 
"Tr0" terms yet avoid double-valuedness on the corners, use a limit procedure in which the 
condition on eg is relaxed in a small neighborhood of the corners. Next, consider the limit 
as this neighborhood "shrinks" to the corners. 

The precise procedure is as follows. Suppose that 65 does indeed coincide with u on t' 
and t", but that on T the tetrad has the form 

Cq = ipu — wip n 

e^ = Tpn — wipu ^ (D2) 

where if) = {1 — w^)^^/^. For each 5 G [0,1], w = w{x;S) is a suitably continuous and 
differentiable point-dependent boost velocity defined on T. Further, for each 6 assume that 



33 



w{x; S) = except on a "small" neighborhood Afs of the corners B' and B". For each S the 
set J\fs is not connected, but is comprised of the disjoint union of two pieces A/5 and Af^'. 
The set Afg is a "small" region of T which contains B', and in the limit 5 — we have that 
(A/] — B') — s> 0. Similarly, the set A/"/' is a "small" region of T which contains B", and in the 
limit 6 ^ we have that (A^" — B") —>■ 0. Finally, for each 6 demand that w{x; 6) = v{x) 
whenever x E B'\J B" . This ensures that on the corner two-surfaces = u and eg = n. 
Our construction provides us with a family of tetrads parameterized by 6. By construction 
the member tetrad corresponding to each value of 6 is TR-gauge on most of T, however, as 
the corners are approached, each member is continuously boosted until it is RT-gauge on 
the corners (see Figure (^). Hence, each 6 tetrad is single- valued on the corners. The idea 
is to use a 6 tetrad in our divergence expression and consider 



-/ d^xJ^TrK = \im-[ d^x 

K JdM <5-»0 K Jm 



pv 



(D3) 



where the expression on the left-hand side symbolically represents the integral of the trace 
of the extrinsic curvature of dM. as embedded in M. over all of dM. (which picks up finite 
corner contributions, since the normal of dM. changes discontinuously from m to on these 
two-surfaces). We can use Stokes' theorem to find 



- / d^x 



-9 



V, ( 



- r d^xVhK + - f d^ 
K Jt' K Jr 



pa 



pX 



(D4) 



Focus attention on the T boundary term 
1 



K Jr 



pu 



K Jr 



(Tx 4-1 np e^^' 



ct' Cp 



p' pa' 



Cp 



(D5) 



We have used the inhomogeneous transformation rule for connection coefficients to express 
the 6 connection coefficients in terms of the connection coefficients F*^ f/,^ determined by the 
TR-gauge tetrad. Using the the boost relations (P32|), we find after some algebra that 



K Jr 



d-^x 



e"" ea " F" 



pu 



1 



K, Jr 



d-^x 



(D6) 



where 9 = — F*^' i-'s-' = — T'^i-'s (s runs over (_L',i,2)) and (f = (p{x;5) = tanh~^(w(x; 5)). 
Next, since u = l/N{d/dt — V), with some integrations by parts the final integral on the 
right-hand side becomes 



1 



— d X ^/^u[ip] = 
K, Jr n Jb 



1 f^" 1 r • \ r / - \ 

d^x^/a(p+- |^d^x^p^/a |^d^x^/aip {6bV''j . 



K Jr 



K Jr 



(D7) 



We have that lim^^o V^(2;; ^) = everywhere on T except for corner points where 
lims^Q(p{x;S) = 0(x). Therefore, in this limit only the first corner-term integrals on the 
right-hand side survive. Hence we have the main result 
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I r I pt" If _ 1 fB" 

K JdM K Jt' K Jr K Jb' 



(D8) 



which justifies Since the action in (|3.5|) is essentially a metric action, we have 

borrowed the results from Ref. to obtain the variation ( |3.6| ). However, it is not difficult 
to use the 6 tetrad method to obtain this result. To perform this calculation it helps to 
assume that 6w = 0, or, in other words, the variations of the 6 tetrad and TR-gauge tetrad 
are "locked" together. We note in passing that a straightforward though somewhat lengthy 
calculation shows that variation of the action (|3.1j ) is 



K JM 
1 



K JM 



(D9) 



One must insert the b tetrad into this expression and then take the limit 5^0. 

The pure imaginary boundary term ( |3.12| ) added to the Goldberg action may also be 
expressed as 



-•5" = 77 , 
The variation of this expression is 



(DIO) 



Ik j m 



d^Xv^V.VA(e'^"^^e^.5e% 



(Dll) 



2k jm 

Using the 5 tetrad in each of the above expressions, one can take the lim^^o and verify that 
—S'^ and —bS'^ contribute no corner terms. 



y.K. . M I \ 
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FIGURES 




FIG. 1. Spatially bounded spacetime region M.. 
In this simple diagram one dimension is suppressed. Although the timelike portion T of the 
boundary of M. is depicted as connected, this is not a necessary restriction. 
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FIG. 2. Radial T foliation of A^. 
The diagram depicts a cut-away view of M. which shows the nested timelike sheets of constant 
The portion of the hypersurface Ef^ shown lies between the two-surfaces -Bt»,rt and -Bt^,r"- Both 
and r* are constants; r" is the constant which labels the actual boundary T = 
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T 

FIG. 3. Geometry of an "undamped" S foliation. 
In this diagram the two-dimensions corresponding to the two-surfaces Bt^r are suppressed (but 
only one dimension of the rather symbolic local lightcones is suppressed). The i = slice as well 
as the final t" slice of the temporal S foliation are depicted as heavy lines. The actual boundary T 
is also shown as a heavy line. Notice that the E Eulerian history (along the integral curves of u) of 
a two-surface B generates a timelike "surface" T which emerges from (or crashes into) the actual 
boundary. The S slices, depicted with dotted lines, are orthogonal to the actual boundary, while 
the S slices are orthogonal to the T sheets. When the clamping assumption is made all S slices 
are also S slices, and, therefore, T and T also coincide. When this assumption is made, we drop 
all bars from the formalism. RT-gauge tetrads have their time leg locked to u and their space leg 
locked to n, while TR-gauge tetrads have their time leg locked to u and their space leg locked to 
n. These vectors are drawn in the diagram as arrows extending orthogonally off Btny/ and Bt^yi. 
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FIG. 4. The 5 tetrad. 

In this diagram the two-dimensions corresponding to the two-surfaces Bt^r are suppressed (but 
only one dimension of the rather symbohc local lightcones is suppressed). The 5 tetrad is TR- 
gauge on T except for small neighborhoods of the corners, for instance, the neighborhood M'^ of 

the corner B^'i ^-f = B" . As depicted in these "snapshots," in such a neighborhood the tetrad is 
continuously boosted as the corner is approached and becomes RT-gaugc at the corner. The dotted 
vectors drawn around the lightcones on A^^' are neither RT gauge nor TR gauge. 
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